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Abstract 

We consider the thermal emission rate of dileptons from a QCD plasma 
in the small invariant mass {Q‘^ ~ but large energy ^ T) range. 

We derive an integral equation which resums multiple scatterings to include 
the LPM effect; it is valid at leading order in the coupling. Then we recast 
it as a differential equation and show a simple algorithm for its solution. We 
present results for dilepton rates at phenomenologically interesting energies and 
invariant masses. 
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1 Introduction 


In the last few years progress has been achieved in understanding the production 
of electromagnetic probes in a quark-gluon plasma in equilibrium. In the case of 
a hard photon (which is defined as having an energy larger than the temperature 
T of the plasma) the usual Compton and annihilation processes discussed long ago 
are not the only leading-order mechanisms. It turns out that bremsstrahlung 
of a quark in a plasma gives a leading contribution and, more importantly, a 
crossed process from bremsstrahlung actually dominates at a large enough energy 
This process, sometimes referred to as off-shell annihilation or annihilation 
with scattering, is of type 3-body —> 2-body : the photon is produced via quark- 
antiquark annihilation where one of the incoming quarks is put off-shell by scattering 
in the plasma. However, the formation time of the photon is comparable to, or larger 
than, the average time between soft scatterings of a quark in the plasma (inverse of 
the damping rate). This implies that multiple scattering plays an important role in 
these photon emission processes 0. 

A correct treatment of multiple scattering, even after performing Hard Thermal 
Loop (HTL) resummation 0, involves the resummation of ladder diagrams with 
effective propagators for both fermions and gluons. This resummation leads to an 
integral equation, which has been solved to give the photon emission rate to leading 
order in the strong coupling Qs [IT,^,^. The resummation of higher loop diagrams 
implements the Landau-Pomeranchuk-Migdal |r^ , p!^ , p!6[| effect which has been much 
discussed recently in the context of jet energy loss in a hot medium ||17| , |I^ ,0 ^ 


However, there are differences between the approach followed here for 
photon production and that used for jet quenching: we work consistently in the 
HTL approach, using the actual distribution of moving charges in the plasma and 
dynamical screening of interactions, rather than a model of static, Debye screened 
charges. Note that no magnetic mass needs to be introduced in order to screen 
transverse gluons, since ultrasoft divergences have been shown to cancel thanks to 
cancellations between diagrams of different topologies 
Preliminary hydrodynamical studies 


(see also 24,25,26 


which take into account the 
above mechanisms of photon production, show that the thermal rate (including the 
contribution from the hot hadronic phase) should clearly emerge above prompt photon 
sources (produced in the early stages of the collisions) in a range from 1 GeV — 3 
GeV to 5 GeV at RHIG and LHG. 

The above developments are quite timely in view of the recent WA98 measure¬ 
ments at SPS 1^^ and of the upcoming results at RHIG where the photon spectrum in 
gold-gold collisions is being measured [^|. There is however an uncertainty concern¬ 


ing the possibility of separating the spectrum of direct photons from that of photons 
which are decay products of resonances, in particular tt^’s. In the energy range of 
interest the 7r° signal may be almost two orders of magnitude above the direct photon 
rate and the subtraction of the vr®, y, ■ ■ ■ background from the data is a formidable 
experimental challenge. 
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In the mean time it is interesting to consider an alternative channel which carries 
the same dynamical information as the real photon production rate, but which suffers 
a different, less important, background. This is the case of the production of low mass 
lepton pairs at large momentum. The invariant mass of the pair should be high enough 
to be above the background of Dalitz pairs from the pions, but small enough that the 
rate remains appreciable. Typically one expects the range 150 Mev < < 500 

MeV to be realistic for our purposes is the invariant mass of the lepton pair). 

Higher mass ranges (above the vector meson masses) should also be explored. 

This rate has already been calculated for the Drell-Yan process , 


for the 2 


processes |]36|J^ , and for bremsstrahlung and off-shell annihilation in the approxima¬ 
tion where the quark undergoes a single scattering in the medium |^. The analysis 
of multiple scattering is clearly called for. Compared to the real photon case, the 
results of [Q show increased technical difficulties associated with the much more 
complicated analytic structure of the diagrams: more unitarity cuts contribute and 
the cancellation of divergences between the various cuts is a subtle affair. In accor¬ 
dance with general theorems of perturbation theory |^, the calculated rate turns 
out to be hnite everywhere except at the qq threshold = 4:M^, M^o the quark 
thermal mass) where there remain integrable (square root and logarithm) singular¬ 
ities. However, this threshold is right where the formation time is the longest, and 
rescattering effects are expected to be largest. Therefore the two loop analysis is 
not reliable here. 

The object of this paper is to compute the imaginary part of the retarded current- 
current correlator, at momentum Q which is hard (g° ^ T) but near the light cone 
at leading order in gg. This current-current correlator is related to the 
dilepton production rate, per lepton species, viaQ 


dN 


t+t- 


d'^xd^Q 


__iniH 

127r4Q2(egYT _ 1) 


( 1 ) 


The calculation requires a resummation of diagrams analogous to that for real photon 
production. We hrst review the integral equation already obtained in the case of a 
real photon (Q^ = 0). Then we extend it to off shell photons. Besides minor kinematic 


changes to the treatment of , this also requires inclusion of longitudinally polarized 
virtual photons. It is then shown how to obtain an ordinary differential equation by 
going to impact parameter space, where the equation is solved numerically. The 
result is more complete, and in some ways simpler, than the two loop perturbative 
one. We hnd that multiple scattering effects completely remove the qq threshold 
behavior found in the perturbative treatment; the dilepton spectrum is nonsingular 
and smooth. Finally, we discuss some phenomenological applications. 


^at leading order in and neglecting the lepton mass 
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Figure 1: Typical diagram which must be resummed to determine ImllJ^. The gluon 
lines are soft and HTL-resummed, the solid (quark) lines are hard, approximately 
on-shell, and resummed to include the dominant imaginary part. 



2 LPM corrections for longitudinal photons 


In this section, we leave aside the production of photons by Compton scattering 
{dQ lO) annihilation {qq —> 'jg) which appear at one loop in the perturbative 
expansion of the HTL effective theory p6| , ^ and we consider only the processes of 
type qi 'jqi and qqi —*• yi where i denotes a quark, an antiquark or a gluon. These 
appear at the two-loop level but because of a strong collinear enhancement, powers 
of T^/Ml 1/gl are generated so that they contribute at the leading order in the 
strong coupling. Higher order diagrams with a ladder topology involve the same 
enhancement mechanism [each new rung in the ladder brings a factor of 0{g‘l/g‘iy\ 
and therefore they also contribute at leading order. The summation of these dominant 


diagrams leads to an integral equation which we now discuss [11,12,13 . 


2.1 Reminder: transverse modes 


In 00,0, Arnold, Moore, and Yaffe summed the set of diagrams of the form shown 
in Fig. |l|, for lightlike external momentum Q. The lightlike condition allowed them to 
consider only the two transverse polarizations; in general H^^" has contributions from 
three polarizations (not four, because = 0), but when = 0 the longitudinal 

polarization also vanishes. In terms of the retarded photon polarization tensor, we 
can rewrite the result presented in these papers as follows: 


1 = 1,2 


e^iVc 

Stt 


^ + 00 


dpo [np{ro) 


^f(Po)] 


(Po’^o)^ 


X Re 


(2vr)2 


2Pn ■ f{P±) 




( 2 ) 


with tq = pQ + qQ, Upipo) = l/(exp(po/T) -|- 1) the Fermi-Dirac statistical weight, and 
where the dimensionless function f{p±) obeys the following integral equation: 

/g2i 

j:^C(h)lf(p±)- f{p± + h)]. (3) 

In this integral equation, 6E = go(pi + Tf,^)/(2po^"o) is an energy denominator arising 
from nearly on-shell quark propagators; it can be interpreted as the inverse formation 
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time of the photon. Here = g^C^T'^/A is the thermal mass of a hard quark 
at leading order (or equivalently, in the HTL approximation). C-p = {N^ — l)/2Nc 
is the Casimir in the fundamental representation of the gauge group SU{Nc)] in 
QCD it is 4/3. The second term arises from the gluon exchange lines in Fig. |l| 
(both those shown and those hidden in the fact that the quark line self-energies 
are resummed). The factor C{1±) is the kernel of a collision integral describing an 
elementary collision of the emitter with another particle of the plasma. In terms of 
the transverse and longitudinal HTL self-energies of the exchanged gluon, 

it is given by 000: 


C{h) = 


dll ^ 


2'k5{Iq 




X 


2Im nQ,(L) 


(L2 - Ren„(L))2 + (Imn„(L))2 


Pr{L)Q,Q, 


(4) 


with Qfj^ = (1, q/q) and (L) the transverse or longitudinal projector for a gluon of 
momentum L. In P], it was shown that this collision kernel can in fact be calculated 
exactly thanks to a sum rule satished by the HTL gluon propagator, leading to a very 
simple expression: 


C{L) 


1 1 


(5) 


where = g‘^T‘^[Nc + Np/2]/3 is the Debye mass in an SU{Nc) gauge theory with 
Np fundamental representation flavors. 

Generalizing the above equations to the case of the transverse modes of a massive 
photon is straightforward, as one needs only to replace in 6E by the following 
combination of the quark mass and the photon invariant mass [^|T^ 




^oo ^ Ks = + —Poro ■ 

% 


(7) 


For > 4M^, this quantity can become negative and 6E can vanish. In this case, 
Mgg should be understood as having an inhnitesimal imaginary part ipo^o^/qo, as one 
can readily see from its derivation from the quark energy denominators. 

This is not enough to fully calculate the dilepton production rate; one must also 
include the longitudinal mode of the off shell (virtual) photon. We calculate the 
contribution of this mode in the next subsection. 

^This substitution is sufficient as long as we are considering hard photons that are not too virtual. 
It can be used if the following inequality is satisfied: 

Q^<^qo- (6) 

If this is not the case, then the kinematics is no longer dominated by small angle scatterings and all of 
the above equations must be reconsidered. Note that this inequality implies that one can interchange 
qo and the 3-momentum q at leading order, in all quantities where no cancellation (go ~ q) occurs. 
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2.2 Longitudinal mode 


The contribution of the longitudinal mode can be derived simply by following a pro¬ 
cedure similar to the one employed in ||TT|,p!^. Let us start by writing: 




( 8 ) 


with the following longitudinal polarization vector 


<(Q) 


(q,0,0,qa) 


(9) 


Making use of the Ward identity 
tensor, we obtain 


njQ) 



= 0 satisfied by the photon polarization 

no°(Q). (10) 


Therefore, we see that we need only the 00 component of the photon polarization 
tensor^. This also makes clear why = 0 for real photons. 

Now we must evaluate n°°((5) summing over all diagrams of the type shown in 
Fig. 0. Recall that the retarded/advanced HTL propagator of a hard quark can be 
approximated by: 




P ± ipoe 


( 11 ) 


where P = {po, \/p‘^ + M^p). It is easy to verify that one can replace the Dirac 
matrix f in the numerator of the effective quark propagator by: 


f e{po) u;(e(po) p) «^(e(po) p) , (12) 

where we denote w{p) = i/j^/p u{p) with u{p) the usual massless^ spinor normalized 
according to u{p)u{p) =1/^ (implicitly, po > 0 in this relation) and ujp = 

Then, one can associate the spinors w and w in the numerator of the propagator with 
the vertices, 

l"" = w{.e{po + /o)(p + /)) 7“" w{e{po)p) , (13) 

if L = (/o,0 is the momentum of the gauge boson entering in the vertex. For the 
vertices coupling the quark loop to the soft gluons, we can simply neglect the gluon 

^Note that, since the sum over polarizations (Q) can be replaced by —g^'' thanks to 

Ward identities, it is the opposite of 11^ that must be added to HP in order to obtain the full 11^^. 

To be consistent, if HTL fermion propagators are necessary then HTL corrections to the qq'y 
vertices are also needed |^. However, the complete effect of this vertex is to ensure that the Ward 
identity is satisfied, by rescaling by an 0((?s) amount with respect to Hoo- By using the Ward 
identity to express H^ in terms of Hoo, we have taken this into account, and can compute Hgo at 
leading order only. 

®This property is true because the HTL resummation preserves the chiral invariance of the bare 
theory. 
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momentum, and write 7 '^ ~ 2P^ with = {po,^{po)p)- It is slightly more compli¬ 
cated to deal with the 'jqq vertices since the photon can be hard, but we need only 
the 0 component 7 °, in the conhguration where the quark and the photon are nearly 
collinear. In this limit, we hnd: 7 ° ~ 2y/pr 2^^/\p(^. 

Now that all the spinor structure has been absorbed in the vertices, performing the 
Dirac’s trace is trivial; Tr(l) = 4. In fact, the vertices coupling the quark loop to the 
gluons are the same as in a scalar theory, and the only differences between scalar and 
fermionic particles he in the coupling to the photon and in the statistical weights. 
Therefore, we can mimic]^ Eq. (^) for the retarded self-energy of the longitudinal 
photon in order to obtain 


+ 00 


ImnjQ) = 


e NcQ f [n^(ro) - 


47r q'^ 


Poro 


X Re 


d^PA 


‘^V\Poro\9{p±) , 


(14) 


where g{p±) is a dimensionless scalar function which describes the resummed coupling 
of the quark to the longitudinal photon, and obeys the following integral equation, in 
complete analogy]] with Eq. (^, 


‘^VlPorol = i6Eg{pj^) + T 


(Ph 

W) 


[g{,PA.) - gi.Pi. + h)] ■ (15) 


Combining the result of and the above formulas for the longitudinal mode, 
we can write the equation for the complete polarization tensor of a massive photon: 


imn.ag) 


271 
X Re 


r’+oo 


dpo [ripiro) - n^ipo)] x 

•00 


(16) 


(27r) 


Pl + rj 

"^iporoP 


P± ■ fiP±) + 


^/\Poro\ Q 


jgiPi 


where /(P±) and gip±) respectively obey Eq. (j^) and Eq. (plsl) . Note that this 
expression is valid for one flavor of the emitting quark, with electrical charge e (but 
Np flavors of quarks are taken into account for the scattering centers, as can be seen 
from the expression of the Debye mass). For the emission by families of quarks 


® The substitution rules are the following: there is a — sign since it is the opposite of 11^ that 
enters in A factor e(poro) conies from the last loop (other factors of e(poro) are kept inside the 
function g{p±)). The factor Q^(q^ comes from the relationship between 11^ and IIoo, and the factor 
(Po + ^o)/(Poro) must be removed from Eq. (j^ since it was there to account for the coupling of the 
quark to the transverse modes of the photon. 

^Just replace the “tree-level” vertex 2pj^ (i.e. the starting point of the iteration) by 2ydporo| 
and the function f{p±) by g{p±)- 
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with electrical charges e*, one must substitute: 



/=i 


( 17 ) 


2.3 Born term 


Before solving numerically the above integral equations, it is interesting to show 
that they contain the contribution of the Born term, corresponding to the direct qq 
annihilation into a virtual photon. Indeed, if we solve the equations and (plSf) at 
order 0 in the collision term, we obtain respectively 


and 


Pi ■ /(pi) = -4p: 


f/(Pi) = -4Vboro| 


2 


go pi + + i 


Porp ^ 

<?0 




go pl + M^g + i 


Porp ^ 
QO 


(18) 


(19) 


These quantities yield a nonzero real part, when integrated over p, provided that 
Mgg < 0. This occurs when Q'^ > 4:M^, i.e. above the threshold for the process 
gg —>■ 7 *. Keeping only the real part and substituting these values in Eq. ( 0 ). we get 


Imn4(Q) = 


e^iVc 

2tt 


b+oo 


dpo^(-^efr) [np{ro) - n^ipo)] 


X 


roPo 


( 20 ) 


which is identical to the expression obtained in for the Born contribution to the 
virtual photon retarded self-energy. Note that the term proportional to can be 
seen as an effect of taking into account the correct HTL vertex correction in the 
calculation in this Born term. In the present approach, we did not need to use these 
vertices explicitly; rather, we relied on them to ensure the Ward identity, used in the 
calculation of the longitudinal contribution (see the derivation of Eq. ([T0|)). 

Therefore, the only leading order terms that are not included in Eq. (pTl) are the 
2 —2 processes considered in This extra contribution should be added to 

our results in order to obtain the complete leading order dilepton spectrum. 


3 Formulation in impact parameter space 

It is possible to approach the integral equations analytically either in an expansion in 
the collision term (the zero and hrst orders are known), or to leading logarithmic order 
in small 6E/C-, but each method is accurate only in a limited kinematic regime. In 
general it is necessary to treat the integral equation numerically. This has been done 














for on-shell photon production, in momentum space, by using a variational approach 

00 - 

However, a simpler method can be obtained after going to impact parameter space 
by a Fourier transform. Doing so is made practical by the fact that the collision term 
is known in closed form and has a rather simple Fourier transform. The method was 
outlined in [§ in the case of real photon production, and can be extended to the case 
of dileptons. 

Let us hrst introduce the Fourier transforms: 

f(pj = j , 

g{pj = f , (21) 


where we use the same letter to denote both a function and its Fourier transform, as 
the argument should make obvious which one enters in a particular expression. The 
functions f{b) and g{b) obey simple partial differential equations: 

i;T^(A4»-Ai)/(6)+gJC,rD(m„6)/(i>) = -2!Vi^(i.) . 

- Ax)g((>) + r/.C,TD(m^b)g{b) = 2^^\S(b) , (22) 

where we denote 


D{m^b) 


1 


7b + ln 



+ Ko{m^b) , 


(23) 


with Kq the modihed Bessel function of the second kind and 7 ^, the Euler constant, 
7 j, ~ 0.577. As we shall see later, the inhomogeneous terms proportional to 5(6) in 
these equations determine uniquely the normalization of the solution by imposing the 
small b behavior of the solution. The quantities that are needed in Eq. (^) in order 
to compute the photon polarization tensor are then given by 

j J^P± ■ fiP±) = Im ■ /( 6 ) , 

Re J g{p^) = £m Re g{b) . (24) 

We must now determine the boundary conditions to impose on the solutions of 
Eqs. (1^) in order to determine them uniquely. First, the functions /(6) and g{b) 
must remain bounded when | 6 | —>■ -|-cx). The most general large b behavior is the 
sum of an exponentially growing and and exponentially shrinking solution; only the 
shrinking solution is allowed, so the required boundary condition is 


lim /( 6 ) = 0 , 

0 —>+00 

lim g{b) = 0 . (25) 

6—>+CXD 


9 














Another boundary condition can be derived at 6 = 0, from the small distance behavior 
of the differential equations satished by f{b) and g{b). It is easy to check that D{m^b) 
vanishes like 6^ ln(6) when b approaches zero, so that the term proportional to D{m^b) 
becomes irrelevant in this limit. In fact, the dominant small b behavior is completely 
determined by the Laplacian and the inhomogeneous term. One hnds easily: 


fib) 


b^0+ 


aib) 


b^0+ 


+ 0 ( 6 ). 

71 qo b 

— ^^yipoi-ol ln(m„6) +0(6”). 


(26) 


Note that even if these small b behaviors are singular, they do not affect Eqs. 
since the singular part of f{b) is purely real, and the singular part of g{b) is purely 
imaginary. Since the differential equations (|22|) are second order affine equations, a 
solution is uniquely determined by two boundary conditions. Therefore, Eqs. (|^) , 
(^) together with the boundary conditions (P3|) and (^UD provide a complete refor¬ 
mulation of the original problem. 

These differential equations can be simplihed further by noting that rotational 
invariance in the transverse plane implies that g{b) depends only on the modulus b, 
and that f{b) can be written as f{b) = bh{b). Therefore, we have 


Ax/(h) = b 


h”{b) + 3 


h\b) 


A^gib) = g"{b) + 


9'ih) 


(27) 


In terms of the new function h(6), the hrst of Eqs. (^) becomes 

(Ppa 


Re 




and the desired differential equation is 


(vi - dl - 3j ) 6(6) + iC,TD(rn^b)h(h) = 0, 


with boundary condition 


Kb) 


2po’^o 

_ 

6^0+ 


+ 0 ( 1 ). 


(28) 


(29) 


(30) 


4 Numerical resolution 

General nonlinear ordinary differential equations with mixed boundary conditions are 
usually solved with the “overshoot-undershoot” method. However, because our differ¬ 
ential equations are linear, they can be solved by a single application of a quadratures 
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ODE solver (specific algorithms can be found in |Q). We will describe the procedure 
for h{b); the procedure for g{b) is completely analogous^. 

What we want is the imaginary part, at the origin, of the exponentially decaying 
solution of Eq. normalized to satisfy the boundary condition at zero, Eq. (|30|) . 
This can be achieved by Ending the decaying solution, with arbitrary normalization, 
and scaling it to have the desired small b behavior (including rotation by a complex 
phase to make the 1/6^ piece purely real). To find the pure decaying solution, it is 
sufficient to begin at large b with arbitrary (nonzero) initial data for h{b) and h'{b) and 
to evolve the differential equation inwards towards 6 = 0. The solution which grows 
exponentially at large b shrinks as we evolve towards the origin, the normalizable 
solution grows. 

To determine what value of b is sufficient to suppress the wrong solution by a 
suitably large factor, we can determine the approximate large b behavior of Eq. (p^ 
by approximating D to be constant (at large b it only varies logarithmically, D{m^b) ~ 
( 7 g + \og{m^b/2))) and dropping the 39^/6 term. The result is, 

h(6>l/m„) ~ c+e^+”^D^ + c_e^-™D\ 

A± = . (31) 

V % 


We should choose initial data at a value of 6, large enough that the ratio of exponen¬ 
tials is abundantly larger than the desired accuracy. It is also easy to 

use the above equations to ensure the initial data for h and h' are almost those of 
the shrinking solution, so the coefficient c+ is initialized close to zero. (This proce¬ 
dure is not exact because of the 39^/6 term we dropped, and the b dependence of D 
which we neglected.) Starting with this initial data, the ODE, Eq. (^), is evolved 
towards 6 = 0 by standard methods, and the coefficients of the 6“^ and 6° behavior 
are extracted. 

For the extraction of the b~^ and 6° behavior, it is helpful to determine the small 
b behavior a little more precisely. At small b, D{m^b) ~ b‘^ln{l/m^b) is negligible; 
but the Mgg term may be large enough, only to be negligible at quite small b, where 
numerical error in the b~^ term can pollute the 6° term. Hence it is useful to match 
to the behavior neglecting the D{m^b) term but including all other terms; 


h{b <C l/'m.jj) 


r\j 


Cl 


Meffiri(Meff6) 

b 


+ C2 


/l(Meff6) 

Mesb 


(32) 


with Ji and Ki the modified Bessel functions of the first and second kind (for < 0 
one should replace K, I with F, J and Mes with |Meff|). The desired quantity, Eq. (p^), 
is just (iporo/Tiqo) Imc 2 /ci. (By working in terms of C 2 /C 1 , we avoid needing to rescale 

®The procedure described in the present paper is slightly different from the one outlined in ||]. 
The latter has also been checked to work, albeit more slowly since it requires two runs of the ODE 
solver instead of one. 
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the solution to the appropriate normalization.) It is also possible to include the small 
effect of D{m^b) perturbatively. 

The treatment of g{h) is completely analogous; the asymptotic large b analysis 
is the same. The differential equation differs only in that the coefficient of the db/b 
term is 1 instead of 3, so the small b behavior is 

g{b < l/mj ~ CiKo{Mesb) + C2lo{Mesb) , (33) 

and the desired quantity, Eq. (^), is {^porQ^/poro/nqo) Imc 2 /ci. 


5 Phenomenology 

The major problem in the 2-loop evaluation of the dilepton rate was the presence of 
a singularity at the location of the tree-level threshold = 4M^. Basically, this 
divergence is a common feature of any fixed loop order calculation near a phase-space 
boundary. A first point to check in the present resummation is that this problem has 
been cured. This is readily seen in the plot of figure In this plot, ‘single scattering’ 
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Figure 2: The sum of all the multiple scattering diagrams compared to the single 
scattering contribution only. In this plot cig = 0.3 and qo/T = 50. 
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denotes the contribution calculated in [^, i.e. what one would get by keeping only 
the first order in an expansion of the integral of eqs. (P) and (|15|) in powers of the 
collision kernel C{1±). One can see that including all orders in C{1±) makes the result 
continuous at = 4M^. Note that the Born term (i.e. the contribution of in 
the solution of the integral equations) is not included in this plot. 

If one now includes also the Born term, an interesting feature of the result is 
that there is no trace of the threshold at = 4M^. In other words, the step 
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function 6{Q‘^ — 4M^) behavior which is characteristic of the Drell-Yan process qq 
7 * ^ is completely washed out when one corrects this process by in-medium 
rescatterings. This is illustrated in hgure This property is easy to understand 



1 10 100 


q2/t2 


Born term 

rescattering corrections - 

2->2 processes . 

Total - 

Figure 3: All the contributions to Imllj^ up to 0{as). In this plot ctg = 0.3 and 
go/T = 50. 


after the integral equations and (lA), whose solution contains the Born term and 
all the rescattering corrections, have been rewritten as differential equations with 
mixed boundary conditions. The parameter enters only as a coefficient in the 
linear differential equations (P^). Because of the term containing D, the solution is 
always regular and exponentially falling at large b (see Eq. (pID). This ensures that 
it evolves smoothly as changes]^. One can also add that the smoothing of the 
threshold depends on the value of the strong coupling constant gs, as illustrated in 
hgure On this plot, one can see a remnant of the threshold at small coupling, while 
it completely disappears at strong coupling. We also illustrate on this plot a simple 
scaling property of the function qo, ,T); from equations (|I6D , (H) and 

( 0 ) one can show that the dependence on as scales according to 


, (34) 

^An equivalent explanation, in momentum space, for the smoothness at threshold is the following: 
in eqs. (i and respectively, the f{p±) and g{pj_) terms under the integral originate from the 
damping rate of the fermion (due to the rescattering of the fermion in the medium) and can be seen 
to provide an imaginary part to 5E. The pole occuring at in the Born term is then 

shifted away from the real axis by an amount proportional to the damping rate and no 0(— 
appears. 
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Figure 4: Imll^ as a function of the photon invariant mass for two values of the 
strong coupling constant. The dotted lines are the corresponding Born terms. The 
circles illustrate the scaling law of Eq. (|3^) . The value of q^/T is set to 30 in this 
plot. 


where F is a dimensionless function of qo/T and (in other words, gT is the 

natural unit for the photon invariant mass in this problem). 

Concerning the respective size of the transverse and longitudinal photon polar¬ 
izations, we have checked that at low Q^/q^, the rate is, as expected, dominated by 
the transverse mode and that the longitudinal mode becomes important only when 
Q'^/Qo approaches 1, a region where rescattering effects are not very large. 

In order to illustrate how the LPM corrections evaluated in the present paper 
affect the dilepton production rate, we present now a series of plots for realistic 
values of the parameters. Figure ^ shows the dependence of the dilepton rate with 
respect to the mass of the lepton pair, the pair total energy being set to 5 GeV, for 
a temperature of 1 GeV (we consider 2 quark flavors, and take Os = 0.3). One can 
see on this plot that considering only the Born term and the 2 —>■ 2 processes is not a 
good approximation. In the low mass region (Q < 1 GeV), the rescattering diagrams 
more than double the contribution of the 2 —> 2 processes. In the intermediate mass 
region, near the tree-level threshold, the effect of these processes is even larger and 
they completely wash out the threshold, as explained before. It is only for large pair 
masses (Q > 3 GeV) that these corrections start being small compared to the Born 
term. However, the ratio Q/qo of the mass over the energy of the pair becomes large 
and the approximations used cease to be reliable. 

One can also study how this mass dependence varies with the temperature. In 
hgure we display the total dilepton rate for two different temperatures (T = 1 GeV 
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rescattering corrections - 

2->2 processes . 
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Figure 5: Dependence of the dilepton rate on the mass of the pair. In this plot 
T = 1 GeV, go = 5 GeV, = 2 and cts = 0.3. 

and T = 500 MeV), the other parameters being the same as in the previous hgure. 
In this plot, the vertical lines indicate how much of the total rate is due to the 
rescattering contributions. In addition to the drop one expects when one lowers the 
temperature, it is possible to see that the effect of these processes are limited to 
smaller values of the mass Q. This is because the scale which controls these effects is 
controlled by thermal masses proportional to T. 

Finally, in hgure |^, we display the dilepton rate as a function of the pair energy, 
the mass Q being hxed. At large energy, one can see the typical exponential drop 
in exp(—go/F). More precisely, the asymptotic behavior of the dilepton rate at large 
go and hxed Q is in y/^exp(—go/T), the factor being characteristic of the LPM 
ehect at large photon energies. In addition, increasing the pair mass has the expected 
ehect of making the rate decrease. Note that the curves in this plot are unreliable 
at the low end, near go ~ Q, since the approximation go ~ g is certainly invalid in 
this region. This is the region where the lepton pair is produced almost at rest in the 
plasma frame (g ~ 0), which precludes any use of a collinear kinematics. 

6 Conclusions 

In this paper we completed the 0{as) calculation of the rate of low-mass lepton pairs 
produced at large momentum in a plasma in equilibrium. Production of the pairs in 
multiple scattering processes are found to be important and even dominate the rate in 
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Figure 6: Dependence of the dilepton rate on the temperature. The vertical lines 
represent the contribution of the processes involving scatterings in the medium. In 
this plot go = 5 GeV, = 2 and ctg = 0.3. 
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Photon energy (GeV) 

Figure 7: Dilepton rate as a fnnction of the pair energy, for Q = 150 MeV, Q = 
500 MeV and Q = 1.5 GeV. In this plot T = 1 GeV, = 2 and cts = 0.3. 
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a wide range of masses. An interesting resnlt is the absence of a threshold in the mass 
dependence of the virtnal photon: the rescattering processes completely wash ont the 
qq annihilation threshold naively expected. The absence of a threshold in the photon 
invariant mass looks different from the recent resnlts of lattice calcnlations 
where a prononnced threshold is observed in the prodnction of a lepton pair at rest in 
qnark-glnon plasma in eqnilibrinm (in the qnenched approximation). We emphasize, 
however, that we cannot directly compare onr resnlts to these lattice resnlts, since 
we consider only low invariant mass dileptons at high energy, that is, close to the 
light cone (in the plasma frame), while the lattice data presented in P3| , |i^ are for 
dileptons at rest-a regime where onr treatment is not valid. 

The present resnlts assnme that the plasma is in eqnilibrinm, and nse asymptotic 
valnes for the various thermal and screening masses. They are therefore expected to 
be quantitatively reliable only at extremely high temperature. More realistic results 
for comparison with present and future experiments at RHIC and LHC could be 
obtained by using more realistic mass estimates, from the lattice for example, and also 
by relaxing the assumption of chemical equilibrium. Of course, the obtained results 
would be on a less hrm theoretical ground but they might be phenomenologically 
useful. 

The dilepton channel can offer complementary information to the real photon 
channel as a probe for formation of a quark-gluon plasma in heavy ion collisions, 
since the background is expected to be relatively much less important than for real 
photon. In this respect it would be useful to evaluate the contribution of non thermal 
sources of lepton pairs for comparison with the rates calculated here. 
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